Neutral-ionic phase transition : a thorough ab-initio study of TTF-CA 
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The prototype compound for the neutral-ionic phase transition, namely TTF-CA, is theoretically 
investigated by first-principles density functional theory calculations. The study is based on three 
neutron diffraction structures collected at 40, 90 and 300 K (Le Cointe et al., Phys. Rev. B 51, 3374 
(1995)). By means of a topological analysis of the total charge densities, we provide a very precise 
picture of intra and inter-chain interactions. Moreover, our calculations reveal that the thermal 
lattice contraction reduces the indirect band gap of this organic semi-conductor in the neutral 
phase, and nearly closes it in the vicinity of the transition temperature. A possible mechanism of 
the neutral-ionic phase transition is discussed. The charge transfer from TTF to CA is also derived 
by using three different techniques. 
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I. INTRODUCTION: 

Charge transfer salts presenting mixed stacks with al- 
ternating donor (D) and acceptor (A) molecules have 
been extensively studied over the last 20 years for their 
original neutral- ionic phase transitions (NIT) i£ Very re- 
cently, this class of materials has gained renewed inter- 
est as it has been demonstrated experimentally that in 
some cases the conversion from the neutral state (N) 
to the ionic state (I) or from I to N can also be in- 
duced by photo- irradiation. 3 Despite intensive theoretical 
work, the mechanism of the phase transition and photo- 
conversion has not yet been clarified. In these systems, 
unlike in other classes of compounds, no dominant in- 
teraction has been evidenced and the nature of the NIT 
must be related to a subtle interplay of different type of 
interactions. 

Depending on the choice of D and A molecules as well 
as on molecular packing, different type of NIT have been 
observed mainly under pressure and sometimes also un- 
der temperature variation. Continuous and discontinu- 
ous NIT have been reported presenting different ampli- 
tudes of charge transfer (CT) variations from D to A, 
often with a dimerization in DA pairs along the stacking 
axis. The TTF-CA complex made from tetrathiafulva- 
lene (D = TTF) and p-chloranil (A = CA) molecules is 
considered as the prototype compound for NIT. At atmo- 
spheric pressure, it undergoes a first order NIT at a crit- 
ical temperature of about 80 K. This symmetry breaking 
phase transition leads to a ferroelectric low temperature 
phase where the initially planar D and A molecules are 
both deformed and displaced to form DA pairs along the 
stacking chains^ According to vibrational spectroscopy 5 
and CT absorption spectral the CT has been estimated 
to be about 0.2 e~ in the neutral high temperature phase 
and 0.7 e~ in the ionic low temperature phase. TTF-CA 
is also a typical example showing N-I and I-N conversion 
under photo-irradiation^ Recently, ultrafast (ps) optical 



switching from I to N has been observed by femtosecond 
reflection spectroscopy^ 

Numerous theoretical approaches have been developed 
to study the NIT. For many years, the balance between 
the energy needed to ionize molecules and the gain in 
the Madelung energy has been at the heart of the con- 
cept of NITAi More advanced models mainly concern 
extensions of the one-dimensional (ID) modified Hub- 
bard model which includes on-site Coulomb repulsion, 
transfer integral and on-site energy. Anusooya-Pati et 
al£ have recently shown that such a model has a con- 
tinuous NIT between a diamagnetic band insulator and 
a paramagnetic Mott insulator. Therefore, the ground 
state at the NIT may be metallic but unconditionally 
unstable to dimerization. Such a metallic behavior near 
the N-I borderline has been indeed reported^ Depend- 
ing on the strength of other interaction such as long 
range Coulomb interaction, electron-phonon (Peierls) or 
electron-molecular- vibration (Holstein) couplings, con- 
tinuous or discontinuous ionicity changes can be gen- 
erated and the dimerization instability is more or less 
affected^ 

Despite increasing computer power, first-principles cal- 
culations are still very rare in the field of molecular crys- 
tals. These calculations offer a unique tool for analyzing, 
at a microscopic level, CT salts without any beforehand 
assumption concerning the relative strengths of electronic 
interactions. They have already contributed to evidence 
the coupling between CT variation and anisotropic three- 
dimensional (3D) lattice contractions in the presence of 
hydrogen bondgifl and to analyze the quantum intra and 
inter-chain interactions in the mixed stack CT crystal of 
TTF-2,5Cl2BQ»iiThis compound has twice as less atoms 
per unit cell as TTF-CA which makes the computational 
effort considerably lower, but only few experimental work 
is available and structural data are limited to the ambient 
condition phase. From an artificial structure obtained by 
relaxing all atomic positions within the experimental unit 
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cell, it has been shown that the weak molecular distor- 
tion due to symmetry breaking induces a non-negligible 
contribution to the total CT variation at the NIT^i 

Starting from first-principles Density Functional The- 
ory (DFT) calculations, we present here a thorough study 
of TTF-CA. Our main purpose is not to describe the 
phase transition itself but to reach a precise understand- 
ing of the electronic ground states on both sides and far 
from the transition point. This work gives a new and 
complementary insight for this family of CT complexes. 
After a brief section concerning computational details, we 
present, in Section lTTTl a complete analysis of the valence 
and conduction bands with the help of a tight-binding 
model fitted to the ab-initio results. This will lead us 
to discuss a possible mechanism of the phase transition 
itself. In Section IIVI the 3D total electron density at 
300 and 40 K will be analyzed in details using Bader's 
topological approach. 12 Section will be devoted to the 
determination of the CT in the high and low tempera- 
ture phases from our tigh-binding model, from Bader's 
approach and from a simple model based on isolated 
molecules calculations. The last method will only be 
briefly mentionned as we defer detailed description to a 
forthcoming publication. 



II. COMPUTATIONAL DETAILS 



We used the projector augmented wave (PAW) 
method 16 which uses augmented plane waves to describe 
the full wave functions and densities without shape ap- 
proximation. The core electrons are described within 
the frozen core approximation. The version of the CP- 
PAW code used for all our calculations considers only 
k points where the wave functions are real. In order 
to increase the number of k points along a*, we had 
to double and treble the unit cell along the a direc- 
tion and our calculations showed that with 3 k points 
along r(0,0, 0) — > X(1,0, 0) results are well converged 
as for TTF-2,5Cl 2 BQii The band structures calcula- 
tions have been carried out with a plane wave cutoff 
of 30 Ry for the wave functions and 120 Ry for the 
densities, using 8 k points: r(0,0,0), X/2(l/2, 0, 0), 
X(l, 0,0), Y(0, 1,0), W/2(l, 1/2,0), W(l, 1,0), Z(0,0,1), 
and S(l, 0, 1) in units of (a*/2, b*/2, c*/2). 

As the topological analysis requires very accurate elec- 
tron densities, we have first verified that results do not 
deteriorate when using only 3 k points along T — > X in- 
stead of 8 and then increased the plane wave cutoff to 
50 Ry for the wave functions. The subsequent topologi- 
cal analysis was performed using the new InteGriTy soft- 
ware package^ which achieves topological analysis fol- 
lowing Bader's approach— on electron densities given on 
3D grids. In order to obtain very accurate results, a grid 
spacing close to 0.10 a.u has been choosen. 



Our ab-initio ground state electronic structure calcu- 
lations (frozen lattice) are based on the experimental 
structures obtained by neutron-scattering experiments at 
300, 90 and 40 Ki4 The symmetry is monoclinic with 
two equivalent chains with alternation of TTF and CA 
molecules along the stacking axis a. This leads to 52 
atoms per unit cell. Above 81 K, the space group is 
P12i/nl and TTF and CA are located on inversion cen- 
ters. At 300 K : a = 7.40, b = 7.62, c = 14.59 A and 
P = 99, 1° and at 90 K : a = 7.22, b = 7.59, c = 14.49 A 
and (3 = 99, 1°. At 40 K, the space group is Plnl with 
a = 7.19, b = 7.54, c = 14.44 A and j3 = 98,6°. TTF 
and CA are slightly distorded and displaced and form DA 
pairs in a ferroelectric arrangement in the a direction. 

We performed electronic structure calculations within 
the framework of DFT using the local density approxi- 
mation (LDA) parametrization by Perdew and Zunger 13 , 
Becke's gradient correction to the exchange energjsi^ and 
Perdew's gradient correction to the correlation energy. 15 
The main difficulty in treating the molecular compounds 
of this family is the presence of non negligible contri- 
butions of dynamical electronic interactions of van der 
Waals type. They are known to be poorly described 
within traditional DFT approximations, more advanced 
treatments being far too time consuming for such com- 
plex systems. Nevertheless, we believe that the presence 
of a significant CT between D and A makes our ab-initio 
results reliable for the determination of the occupied elec- 
tronic states, especially in the I phase where electrostatic 
interactions largely dominate. 



III. AB-INITIO BAND STRUCTURE 

Our ab-initio DFT calculations in the crystal pro- 
vide not only the total electron density n(r) but also 
for each band the energies and wave functions at all k 
points given in section ^ The present section is de- 
voted to the sole valence bands (VB) and conduction 
bands (CB) which are directly related to the strong 
anisotropy observed experimentally in transport and ex- 
citation properties ii2ii2i2£ The detailed analysis of n(r) 
will be given in section llVl 



A. Quasi ID shape of valence and conduction 
bands 

The dispersion curves for the VB and CB are given 
in Fig. n for the experimental structures at 300 K (high 
temperature phase) and 40 K (low temperature phase). 
In both cases, they are separated by more than 1 eV 
from the other occupied and unoccupied bands forming 
thus a nearly isolated four-band system. The presence of 
two equivalent chains in the unit cell related by a gliding 
plane makes the VB and CB twofold, splittings in some 
particular directions being due to small interactions be- 
tween the symmetry related chains. The dispersion is 
maximum along T — > X which is the direction of recipro- 
cal space corresponding to the chain axis ai2i 

The isodensity representation of the VB states are 
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FIG. 1: Valence and conduction bands of TTF-CA in the 
high temperature phase (300 K, on the left) and low tem- 
perature phase (40 K, on the right). Dots correspond to the 
ab-initio values and curves to the tight-binding model pre- 
sented in Section IhTI 



given in Fig. [3 At 300 K for k = T (Fig. Hi), the density 
is located on TTF with a shape very close to that of the 
highest occupied molecular orbital (HOMO) of TTF& 
At k = X, (Fig. 03), hybridization occurs and the den- 
sity centered on the acceptor molecule is similar to the 
one of the lowest unoccupied molecular orbital (LUMO) 
of CA. 2 2, The same happens in TTF-2,5C1 2 BQ 21 and is a 
direct consequence of the presence of an inversion center 
in the crystal combined to different symmetries of the 
molecular orbitals which are involved. As the inversion 
center is lost at 40 K, hybridization is no more symmetry 
forbidden at k = T and the pairing of TTF and CA is 
clearly evidenced in Fig.|2t (DA pair on the bottom). On 
these figures one can identify the regions of wave function 
overlap responsible for the intra-chain hopping integrals. 

The dispersion is much less important along T — > Y, 
nearly vanishing in the third direction and very simi- 
lar for both structures. As for TTF-2,5Cl 2 BQ f ii the 
band gap is most probably indirect: E g — 0.05 eV at 
300 K and 0.15 eV at 40 K. It is well known that the 
LDA used in our calculations underestimate the gap of 
semiconductors; these values are thus lower estimates of 
the true band gap. However, because of missing data 
a direct comparison to experiments is impossible. Opti- 
cal absorption measurements 2 " showed a first broad peak 
around 0.7 eV for the N phase which is attributed to the 
ionization of one DA pair. It may correspond either to 
band to band transitions or to an exciton peak. The op- 
tical spectra of TTF-2,5Cl2BQ 2 shows a shift to higher 
energies which is coherent with a larger calculated band 
gap for TTF-2,5C1 2 BQ: 0.14 eV in the high temperature 
structured 



B. Tight-binding model 

As shown in Fig. |2J the ab-initio VB and CB can be 
interpreted in the frame of a tight-binding model as lin- 
ear combinations of the HOMO of TTF and the LUMO 
of CA. One should emphasize that in such a model, hy- 



bridizations with lower bands are completely neglected. 
The interaction with deeper occupied states is neverthe- 
less indirectly taken into account via the values of the 
parameters, which are fitted to the VB and CB given by 
an ab-initio calculation of all states. We will follow a 
similar procedure to the one used for TTF-2,5Cl2BQ, n 
but here the task is more involved as we have to handle 
four bands instead of two. Let us write |r, D) and |r', A) 
respectively the HOMO of TTF located at position r and 
the LUMO of CA at r'. The VB Bloch functions are then 
defined by: 



l* lk ) = E eIkR " ( C dI R - D) + Cft |R„ + v, D'> 

n ^ 

+ C£|R n + T,A) + C£|R n + T + i/,A')Y (1) 



Here 1 = 1 , 2 is the band index; R n corresponds to 
a primitive translation; r = a/2 represents a transla- 
tion between A and D along the chain axis and v — 
a/2 + b/2 + c/2 a translation between the two symmetry 
equivalent chains, the molecules from the second chain 
being denoted by D' and A'. The molecular orbitals |D), 
|A), |D') and |A') are calculated using ab-initio simu- 
lations for isolated molecules in their experimental crys- 
talline conformations. For each k point, the weight |C- k | 2 
(j = A or D) is obtained by mean square minimization in 
the unit cell of the difference between the ab-initio value 
of |\& lk (r)| 2 calculated in the crystal and the one deduced 
from expression Q). Fig. [3] displays these weights at the 
6 corners of the reduced Brillouin zone used in these cal- 
culations for both the high (left) and low (right) tem- 
perature phases. In the former, the valence band at the 
r is a pure TTF state for symmetry reasons: C^ = 0. 
The band gap separates nearly pure D occupied states 
from nearly pure A empty states in the whole Z — U — Y— r 
plane, whereas electronic states are of mixed D-A nature 
in the X — W — T — S plane. In the latter, a significant 
mixing of D-A states takes place in the whole Brillouin 
zone, for the VB as well as for the CB. 

The charge transfer pvB from TTF to CA is given by 
the mean value of the weights of the acceptors in the VB: 



PVB 



^E(<l C An + <|C^I 

8WzbV 



(2) 



Fig. shows clearly that the increase of pvB in the low 
temperature phase results mainly from additionnal hy- 
bridization, which is no longer symmetry forbidden, in 
the whole T-Y-U-Z plane. 

In order to model VB and CB in terms 
of a tight-binding scheme, we define the on- 
site energies by (D|H|D)= (D'|H|D'}= and 
(A|H|A)= (A'|H|A')= E£. Two families of inter- 
molecular interactions have been defined, the first 



r at 300 K 



X at 300 K 



r at 40 K 



FIG. 2: Isodensity representation of the valence band states for one chain (contours at 0.0004 e"/a.u. 3 ) of TTF-CA: (a) for 
k = r (left) and (b) for k = X (middle) in the high temperature phase (300 K) and (c) for k = V (right) in the low temperature 
phase (40 K). 
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FIG. 3: Weights on the acceptor |Ca| 2 of the occupied va- 
lence states in the high temperature phase (300 K, on the left) 
and low temperature phase (40 K, on the right). Dots corre- 
spond to the ab-initio values and curves to the tight-binding 
model. 



corresponding to interactions between molecules of 
the same species shown on Fig. 0] and the second to 
interaction between TTF and CA shown on Fig. and 
Fig. Due to the different symmetries of the HOMO 
of TTF and LUMO of CA, the sign of the hopping 
integrals between TTF and CA is alternately plus and 
minus. For example, along the stacking chains one has 
(0,D|H|a/2,A) = t and (a/2, A|H|a, D) = -t. The 
symmetry breaking in the low temperature phase leads 
to additionnal terms so that the latter hopping integrals 
become respectively t + e and — t + e whereas ±t;, ±9i 
and (5- (i = 1 or 2 and j = D or A) become ±tj + e- u 
±9i + Ai and 6[ ± v\ respectively. 

Using the ab-initio energies and the weights |Cj k | 2 at 
the 6 corners of the reduced Brillouin zone used in these 
calculations, we determine the parameters given on Ta- 




FIG. 4: Definition of the interaction parameters between two 
molecules (denoted by j, j= D or A) of the same species. 




2b 



FIG. 5: Definition of the interaction parameters between 
TTF and CA in the (a, b) plane. ±t is the hopping integral 
along the 7r chains. The orientation of the molecules in the 
chains parallel to the a/2 + b and —a/2 + b being different, 
the hopping integrals ti and t2 are also different. 



blcH]to|nT]and the curves plotted in Fig. □ and Fig. 
Some parameters like (3[ and j3 3 2 (Fig- El cannot be deter- 
mined separately because of the small number of k points 
available, but all these contributions are found to be neg- 
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TABLE III: Additionnal deformation parameters (eV) due to 
the symmetry breaking in the low temperature phase (40 K) 
as defined in Section Till Bl 
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0.061 
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-0.003 


-0.004 



FIG. 6: Definition of the interaction parameters between 
between TTF and CA in the (b, c) plane. 



TABLE I: On site energies and interaction parameters (eV) 
between two molecules of the same species as defined in Fig. 2] 



300 K 



40 K 



300 K 



40 K 



/3 b ° 
5? 



2.894 
0.002 
-0.022 
0.003 
-0.002 
-0.004 



3.191 
0.001 
-0.005 
0.004 
-0.003 
-0.005 



Pi 

A A + /3 2 A 

5t 



3.048 
-0.017 
-0.023 

0.002 
-0.004 
-0.002 



3.348 
-0.024 
-0.003 

0.003 
-0.004 
-0.003 



ligible. In Fig. ^ the degeneracy is lifted at T due to the 
small inter-chain interactions 5\ (Fig. 0}, at S due to 9\ 
(Fig.[5J) and, in the low temperature phase at Y due to v\. 
The 6*j are also responsible for the small HOMO-LUMO 
hybridization at Y and Z (Fig.[3J). 

The transfer integral ±t along the chain axis (Fig. 
dominates and its estimate is in agreement with previous 



ones 



6,18.21,24 



In the low temperature phase, t increases 
significantly (Table [H| partly due to the lattice contrac- 
tion and partly to the symmetry breakingiAI As for TTF- 
2,5Cl2BQ, the main effect of the latter is expressed by 
the additional parameter e which is as large as a third of 
t (Table HTD and which leads to a large HOMO-LUMO 
hybridization in the Z — U — Y— r plane. The other addi- 
tional parameters Ai and t>- have negligible values. 



C. Towards a metallic state? 

Our calculations at 300 and 40 K present the electronic 
structure far from the transition. In order to understand 
how the band structure is changed by thermal contrac- 
tion, we performed also a calculation just above the NIT 
for the crystal structure known at 90 K. The correspond- 



TABLE II: Interaction parameters (eV) between TTF and 
CA as defined in Fig. |5] and Fig. |nj 





300 K 


40 K 


t 


0.167 


0.206 


tl - 1 2 


0.003 


0.005 


0i 


0.016 


0.012 


82 


0.005 


0.010 



ing dispersion curves are given on Fig.[7| They display a 
closed indirect band gap between T and Y, with nearly 
pure D occupied states at Y at a higher energy than 
nearly pure A empty states at T, as our method imposes a 
fully occupied VB. This is of course unphysical and is due 
to the DFT-LDA error on the band gap value. It leads 
however to interesting conclusions: (i) the thermal lattice 
contraction reduces the indirect band gap in the neutral 
phase, going toward a metallic state and (ii) as this band 
gap closes, an electron transfer occurs necessarily from D 
to A in the vicinity of the transition temperature, leading 
to an increase of the ionicity of the molecules, even in a 
frozen lattice. 



Within our tight-binding model, the band gap decrease 
results from (i) the decrease of the difference — E° be- 
tween the on site energies, (ii) the increase of the inter- 
action between two TTF molecules along b (|/3°|) which 
increases the energy of the valence band at Y and (iii) the 
increase of the interaction between the two CA molecules 
along a |) which decreases the energy of the CB at 
r (Table ITVjl . The cell contraction increases also t re- 
sulting in a slight increase of the weights of the acceptors 
in the VB |C^| 2 in the X-W-T-S plane from 0.36 to 
0.38 and a very small variation of pvb as observed ex- 
perimentally 6 before the transition. One can notice that 
|C^| 2 in the X — W — T — S plane has already reached at 
90 K the value that we obtained at 40 K. 

Of course, in the actual crystal, the charge transfer 
is largely facilitated by dynamical effects in the vicin- 
ity of the transition temperature, when the gap be- 
comes narrow. Thermal ID excitations from D to A 
have already been observed^ They are precursors of the 
global changes in the crystal which certainly involve intra 
and intermolecular vibrations. Calculations on a A-D-A 
model comple>*ifi have recently shown that the CT varia- 
tion occuring in the VB (along the mixed stack 7r-chains) 
induces an intra-molecular electronic redistribution, af- 
fecting deep molecular states as far as 10 eV below the 
frontier orbitals. These deep states are also concerned 
in the hydrogen bonds (along OH-chains resulting from 
coplanar alternating D and A molecules in the b±a/2 di- 
rection) so that their modifications are directly related to 
the lateral lattice contractions. In the low temperature 
phase, the gap opens again between mixed D-A states 
both in the VB and CB, in the whole Brillouin zone. 
This suggest an electronic mechanism for the NIT where 
the electronic redistribution associated to the structural 
phase transition allows the crystal to avoid the metallic 
state. 8 
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FIG. 7: Valence and Conduction bands of TTF-CA for the 
experimental structure at 90 K. Dots correspond to the ab- 
initio values and curves to the tight-binding model presented 
in Section ITTT1 



combined with the local form of the Virial theorem 12 , 
the local contribution to the potential energy den- 
sity V(r C p) = 1/4 x V 2 n(r CP ) - 2G(r CP ) which have al- 
ready been used to study hydrogen bond strengths^ In 
the present work, V(rcp) will be used to quantify the 
strength (intensity) of the intermolecular interactions. 
It allows to quantitatively compare interactions between 
two or more atoms. This procedure is less speculative 
than the comparison of interatomic distances to the sum 
of van der Waals radiii, or related criteri a 28 i 29 , as van der 
Waals radii are not defined with sufficient accuracy^Si 3 ^^ 
and fail to take into account different relative orienta- 
tions. 



TABLE IV: Parameters of the tight-binding model (eV) 
mostly affected by the structural changes occuring between 
300 and 90 K. 





300 K 


90 K 


n, — tj 


0.154 


0.135 




-0.022 


-0.028 




-0.017 


-0.022 


t 


0.167 


0.189 



IV. INTERMOLECULAR INTERACTIONS 

The previous section showed that VB and CB are 
clearly ID. The aim of the present section is to explore 
thoroughly the total electron density n(r) given by our 
first-principle calculations. This is essential to elucidate 
all features, not only the ID aspects, presented by this 
NIT. 

We could plot isodensity representations of n(r) in 
some selected planes to evidence the most important in- 
termolecular couplings. But, the presence of two chains 
per unit cell related by a gliding plane makes it an end- 
less task and Bader's approach^ is a clever manner to 
circumvent this difficulty 



A. Bader's theory 

Within the Quantum Theory of Atoms in Molecules^ 
n(r) can be analysed in details by means of its topo- 
logical properties. The topological features of n(r) 
are characterized by analysing its gradient vector field 
Vn(r). Here we will focus on bond critical points 
(CP) where Vn(r) vanishes. They are characterized 
by the density at the CP n(rcp), a positive curva- 
ture (A3) parallel and two negative curvatures (Ai,Aa) 
perpendicular to the bond path. The ellipticity of a 
bond, denned as e = A1/A2 — 1, describes the devi- 
ation from cylindrically symmetric bonds. For closed 
shell interactions, further information is obtained by us- 
ing Abramov's 2 ^ expression for the kinetic energy density 
G(r C p) = 3/10(37r 2 ) 2 / 3 n(r C p) 5 / 3 + l/6V 2 n(r CP ), and, 



B. Intra-chain interactions 

In the high temperature phase (300 K), the strongest 
contact within the chains is the one connecting S4 to 
C13-C15 (Fig. |S]and Table E) . It corresponds to the one 
already visible in the VB (Fig. and has thus a sig- 
nificant contribution resulting from HOMO-LUMO over- 
lap. It has a quite large ellipticity which is related to the 
fact that interaction do not occur between two atoms but 
rather between S4 and the bridge of C13-C15. Such a 2D 
attractor is also found for the interaction Ci6 • • ■ C1-C2 
which is 40 % less intense than the former and somewhat 
less than Si • • ■ Cl 2 , Si ■ • ■ Cl 4 , Cl 3 • • • C 4 and Cli • • • C 6 . 
These four latter atom-atom interactions show that CI 
plays an important role for the electronic coupling along 
the chains. The acceptor in the middle (left) of Fig. [S] 
is coupled to TTF on bottom (top) by Cb, Cl 3 • • • C 6 , C 4 
(S 2 ) on the left side and Cl 2 ,Cl 4 ---Si (C 5 ,C 3 ) on the 
right side. 

In the low temperature phase (40 K), the pairing of 
the molecules in DA pairs is clearly visible on these lo- 
cal contacts (bottom of the figure on the right part of 
Fig. |S]and Table Ivjl . The S4 ■ ■ • C13-C15 contact remains 
the dominant one with a 50 % increase of its strength. 
Due to molecular reorientation, the contact Ci6 ■ • • C1-C2 
has changed into a more atom-atom like contact between 
C 2 and C14. Its intensity has more than doubled and 
with S4 • • • C13-C15 they dominate all other local con- 
tacts. The latter have also increased by about 35 % 
and additionnal strong contacts 2 • ■ ■ S 2 and S3 ■ • • Cl 2 
have appeared. These contacts confirm- 1 that inside the 
pair, the CA molecule is essentially connected to the left 
part of the TTF molecule below. Between the DA pairs 
(top of the right part of Fig. |SJ, we observe a small de- 
crease and a more atom-atom like contact for S3 • • • C14- 
Ci6 whereas the S • • • CI and CI • • • C contacts moderately 
increase, CI4 ■ • ■ C3 being the strongest one. 



C. Inter-chain interactions 

The inter-chain interactions can be split in two groups, 
the first corresponding to interactions between two chains 
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300 K 40 K 



FIG. 8: Intra-chain interactions in the high temperature phase (300 K, on the left) and low temperature phase (40 K, on the 
right). Diamonds correspond to the bond critical points and lines indicate the bond path associated to each critical point. 

tions are diatomic. On Fig. EI it is clear that CI4 interacts 
with the C - H 2 bond. In TTF-CA, the importance of 
hydrogen bonds has already been underlined 4ii£L2iL Our 
results show clearly that they dominate all other intra 
and inter-chain interactions in both high and low tem- 
perature phases. Oi ■ ■ • H3 in the OH chains is stronger 
than Oi • • • H2 as it is shorter and more straight, the CA 
and TTF molecular planes being aligned in the b ± a/2 
direction. The other strong interactions reported in Ta- 
ble have comparable V(rcp)'s to the intra-chain ones 
given in Table |V| but there are fewer of them in the same 
energy range. 

In the low temperature phase, the loss of inversion 
symmetry leads to two inequivalent distances for each 
contact. As already discussed by Le Cointe et ali^, some 
of the distance reductions (Table IVI|) are larger than 
those expected from the thermal unit cell contraction 
which is about 1 % along b and c and 3 % along a di- 
rections. The hydrogen bonds are the most affected and 
become 80 % and 50 % stronger for the short and long 
bonds respectively. All other interactions except the long 
S • • • H2 bond, undergo a smaller strengthening. 



D. Discussion 

Most of the contacts discussed in the Section llV Bl and 
from each other to allow direct atom-atom interactions. IIV CI have already been pointed out on the basis of a 
In Fig. E|and Table IVT1 we have reported the properties of detailed comparison between interactomic distances and 
the strongest inter-chain contacts. For all of them except van der Waals radii«^> Our results demonstrate that a 
the two CI ■ • • H contacts, e is quite small and the interac- topological analysis of the total charge density allows 



TABLE V: Strongest intra-chain contacts for high (300 K) 
and low temperature (40 K) phases, e and V(rcp) give respec- 
tively the ellipticity and potential energy density (kj. mol~ ) 
at each bond critical point. At 40 K, intra and inter stend 
for interactions inside a DA pair and between two DA pairs 
respectively. 











300 K 


40 K intra 


40 K inter 








e 


V(rcp) 


e 


V(rcp) 


e 


V(rcp) 


s 4 -- 


•C13- 


C15 


7.4 


-9.3 


5.8 


-16.5 






S 3 -- 


■ C14- 


Cl6 


7.4 


-9.3 






4.4 


-7.8 


Si-- 


Cl 2 




1.3 


-7.0 


1.6 


-10.6 






s 2 -- 


•Clj 




1.3 


-7.0 






1.3 


-7.8 


Cla- 


••C 4 




1.3 


-6.7 


1.3 


-8.9 






Cl 4 - 


• c 3 




1.3 


-6.7 






1.2 


-8.7 


o 2 - 


•Sa 








2.1 


-8.9 






Cli- 


••Co 




2.2 


-6.0 


2.5 


-8.1 






Cl a - 


••c 5 




2.2 


-6.0 






1.9 


-6.8 


Ss-- 


•Cla 








2.2 


-8.1 






Si-- 


■Ck 




1.3 


-6.0 


1.6 


-7.8 






S 2 -- 


Cl 3 




1.3 


-6.0 






1.1 


-7.1 


C16 • 


••Ci 


-Ca 


6.0 


-5.6 










C14 • 


••C 2 








4.5 


-11.7 







separated by a primitive translation along b and the sec- 
ond to the interactions between two chains related by the 
gliding plane. In the c direction the chains are too far 
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x = 







\ X CA=- 1/2 ^ 




X JTF = 






v\ 4 T 



H. 



X TTF=° 



x JTF =-1/2 




FIG. 9: Inter-chain interactions in the high temperature phase (300 K). Diamons correspond to the bond critical points and 
lines indicate the bond path associated to each critical point. 



TABLE VI: Strongest inter-chain contacts for high (300 K) 
and low temperature (40 K) phases. V(rcp) corresponds to 
the potential energy density (kj. mol -1 ) at the bond critical 
point, d is the distance in A between the two atoms of the 
first column. Ad gives the relative variation of this distance 
with respect to its value at 300 K. 







300 K 


40 K short 


40 K long 






V(rcp) 


d 


V(rcp) 


Ad 


V(rcp) Ad 


Oi- 


•H 3 


-16.3 


2.35 


-26.7 


-7% 


-24.3 -5 % 


Ss-- 


•Ch 


-6.8 


3.57 


-10.3 


-5 % 


-8.7 -3 % 


Cli- 


•H 3 


-6.7 


3.03 


-10.8 


-6 % 


-8.0 -1 % 


Oi- 


•H 2 


-12.0 


2.50 


-21.4 


-8 % 


-18.3 -6 % 


cv 


■ ■ C4 — H2 


-8.8 


3.53 


-11.6 


-3 % 


-11.0 -2 % 


Si-- 


■H 2 


-7.8 


2.93 


-10.4 


-4 % 


-7.1 +1 % 


Cl 2 - 


■•Cl 3 


-7.5 


3.53 


-9.1 


-2 % 


-8.1 -1 % 



to go a step further. It leads to a very precise picture 
of the atoms involved in the different contacts (bond 
path) and provides a quantative way to compare their 
shape (e) and intensity (V(rcp)). From Table |VT1 for 
the inter-chain interactions and Table and Ref. |4J for 
the intra-chain ones, we find the following relation be- 
tween the relative variations of the potential energy den- 
sity AV(r C p) = (V40K - V 3 ook)/V 3 ook and the length 
of the bonds Ad = (d 40K - d 30 oK)/d 30 oK: 



AV(r CP ) ~ -10 x Ad. 



(3) 



Thus, the strengthening (diminution) of a given contact 
defined by V(rcp) is directly correlated to the reduction 
(augmentation) in distance. 

All our results on TTF-CA show that HOMO-LUMO 
overlap and hydrogen bonds dominate the intermolecular 
interactions in both high and low temperature structures 
and drive the molecular deformation and reorientation 
occuring at the NIT. The other inter-chain contacts do 
not show any peculiar behavior and are dragged along 
with the former. In addition, many strong interactions 
inside and between the chains involve CI atoms which 
must largely contribute to stabilize these crystal struc- 
tures. 



V. CHARGE TRANSFER 

CT is a key factor of the NIT but it is is well known 
that it is not a well defined quantity and there are numer- 
ous ways to extract it from experimental or theoretical 
results. In TTF-CA it has already been estimated from 
the intensity of the CT absorption band 6 , from the de- 
pendence on molecular ionicity of either bond lengths^ or 
vibrational frequencies^ In this section, we will use three 
different techniques to estimate the CT between TTF and 
CA. The first two are based on our ab-initio calculations 
in the crystal while the third is obtained from a model 



9 



TABLE VII: Estimates of the charge transfer (CT) from TTF 
to CA. 



Tight 


Tight 


Integration 


Isolated 


exp. a 


Binding 


Binding 


Atomic 


molecules 




E g 


E g + 0.5 eV 


Basin 


Model 




300 K 0.59 


0.43 


0.48 


0.47 


0.2±0.1 


40 K 0.72 


0.69 


0.64 


0.60 


0.7±0.1 



Prom Ref. 



build on isolated molecules calculations. 



A. From the tight-binding model 

Within our tight-binding model, the CT in the VB is 
obtained by integrating the analytic expression of IC^I 2 
(Eq. (J2J). We have first estimated pvb by using the 
parameters determined in Section IIII Bl leaving out all 
those which are smaller than 0.005 eV in both phases. A 
second estimate has been obtained by keeping only the 
parameters related to the dispersion along T — ► X. In 
that case, |C^ | 2 is given by 

ICaT = \ (4) 

(E A - E D ) 

2^(E A - E D ) 2 + 4t 2 sin 2 ^f + 4e 2 cos 2 ^ 

Both approaches lead to the same values of pvb indicated 
in Table IVlTl where the corresponding experimental val- 
ues have also been reported in the last column. Even 
though experiments lead to coarse estimates of the CT 
as its determination is always indirect, our high temper- 
ature value seems too large. Two reasons can be put 
forward to explain this discrepancy: (i) the tight-binding 
model is based on a linear combination of molecular or- 
bitals of isolated molecules which disregards their defor- 
mation due to intcrmolccular interactions in the crystal 
and (ii) the underestimate of the dynamical part of the 
electron-electron interactions due to the LDA approxima- 
tion. Among other things, the later will affect the band 
gap which is directly included via E D - E A in Eq. 0] In 
order to get a feeling for the influence of the band gap on 
Pvb, we increased it by 0.5 eV in both phases to locate it 
approximately near the first absorption peak^ The cor- 
responding values have also been reported in Table IVTTI 
In the high temperature phase, it produces a significant 
decrease of 0.16 e~ whereas it is only of 0.03 e~ in the 
low temperature phase. Thus pvb as well as its variation 
between 300 and 40 K becomes in much better agreement 
with the experimental results. 



B. From Bader's theory 

Within Bader's approach, a basin can be uniquely as- 
sociated to each atom. It is defined as the region con- 
taining all gradient paths terminating at the atom. The 
boundaries of this basin are never crossed by any gradi- 
ent vector trajectory and atomic moments are obtained 
by integration over the whole basinpi2*i£ The CT between 
TTF and CA is obtained by summing the atomic charges 
belonging to each molecule. For this purpose very high 
precision is needed as we are looking for a few tenths of 
an electron out of a total of 116 valence electrons. The 
corresponding CT are reported in Table IVTTI These es- 
timates are no more directly affected by the band gap 
error due to the LDA as they are based on the sole occu- 
pied states. Accordingly, they are both at 300 and 40 K 
smaller than those obtained by the tight-binding model 
and thus closer to the experimental results. Neverthe- 
less, in the high temperature phase, the CT seems still 
somewhat overestimated. 



C. From a molecular model 

We obtained a third estimate of the CT by means of 
a simple model based on isolated molecules calculations 
which will only be briefly reported in the present pa- 
per. Here, we are looking for the charge transfer p which 
minimizes the total energy in a fixed structure, i.e. ei- 
ther in the high or low temperature structure. The first 
contribution to the variation of the total energy is the 
molecular one (Carloni et ali^i): 

AE mol (p) = p(e° A - e° ) + ip 2 (U D + U A ), (5) 

where the reference state is the one without any CT, 
Cj is the energy of the relevant orbital in the reference 
state and Uj the Coulomb repulsion of an electron in 
that orbital. The second contribution is the electrostatic 
one limited to the Madelung contribution qiQj/rij. 
The atomic charges qi and qj are deduced from PAW 
calculations on isolated molecules by using a model 
density reproducing the multipole moments of the true 
molecular charge density*2si The remaining contributions 
to the total energy are usually accounted for by a van 
der Waals term which is constant in a given structure 
and which can thus be disregarded. The CT values given 
by this simple model are also given in Table IVTTI 



D. Discussion 

It is satisfying that all our three estimates are in 
close agreement with each other in both low and high 
temperature phases. At 40 K, they are also in good 
agreement with the experimental value. At 300 K, 
we observe a systematic overestimate of about 0.2 e~ 
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with respect to the values obtained by Jacobsen et 
al£, our result being closer to the one deduced from 
intermolecular bond lengths (0.4 e - )^ 

This discrepancy may have different origins. First of 
all, the CT is not a uniquely defined quantity as its eval- 
uation relies on different (and not necessarily equivalent) 
models based cither on theoretical or on experimental 
results. Starting from experimental data, the models can 
only be very simple and most of the results published 
in the literature are based on the behavior of isolated 
molecules upon ionization: bond-lengths, intramolecular 
vibration frequencies, optical spectra..., neglecting thus 
all crystalline effects apart from the charge transfer on 
these molecular properties. 

On the other hand, our theoretical results are based on 
the determination of the charge density. They are all af- 
fected by the error due to the LDA approximation to the 
exchange-correlation potential. It is well known that this 
approximation, even in its gradient-corrected version, un- 
derestimates the dynamical part of the electron-electron 
interactions. In weakly interacting systems where these 
dynamical effects (or so-called van der Waals interac- 
tions) are dominant, this leads to an overestimate of 
the electron derealization and an overbinding effect, to- 
gether with an underestimate of the band-gap. An over- 
estimate of the intermolecular CT is thus expected. How- 
ever, it has already been showniS that in the ionic phase 
the multipolar static effects are considerably stronger 
than in the neutral phase, so that the LDA error on the 
CT should also be considerably reduced in the low tem- 
perature phase. 

In our tight-binding model we have clearly evidenced 
the effect of the underestimate of the band gap. Within 
Bader's approach, the CT determination only relies on 
the occupied states, but leads nevertheless to an overes- 
timate in the neutral phase. Regarding our model based 
on isolated molecules calculations, it has been shown that 
the ionization energy of TTF— and the electron affin- 
ity of CAS 3 - are correctly obtained within the gradient- 
corrected LDA approximation used in the present study. 
But our model takes into account neither the atomic mul- 
tipolar contributions nor the polarization due to the crys- 
tal field. 

Thus, our three estimates are coherent and lead prob- 
ably to a correct determination of the CT in the ionic 
phase, and a too large one in the neutral phase. At 
this stage however, the discrepancy with the available 
experimental values has not to be attributed to the sole 
LDA problem, but also to the definition itself of this CT. 
A good way to check the rather high value obtained in 
this work at 300 K is to perform an experimental charge 
density acquisition by X-ray diffraction and use Bader's 
approach to extract an experimental CT value coherent 
with our theoretical estimate. 



VI. CONCLUSION 

This is the first extensive ab-initio study of the elec- 
tronic ground states of TTF-CA, prototype compound 
of the NIT. We should remind here that temperature is 
included only via the experimental structures, all calcu- 
lations being performed in a frozen lattice. 

The valence and conduction bands are dominantly of 
ID character along the stacking chains and are essen- 
tially a linear combination of the HOMO of D and the 
LUMO of A as in TTF-2,5C1 2 BQ M. The weak molecular 
distorsion due to the symmetry breaking part of the NIT 
is responsible for the main contribution to the charge 
transfer variation. We found also that the thermal lat- 
tice contraction leads to a decrease of the band gap in 
the neutral phase, going toward a metallic state. This is 
consistent with the theoretical work of Anusooya-Pati et 
al£ who found a metallic ground state at the transition 
which is unstable to dimerization. It has also to be re- 
lated with the metallic behavior observed by Saito et al£ 
in an other DA alternating CT complex. 

A detailed analysis of all intermolecular bond critical 
points shows that molecular deformation and reorienta- 
tion occuring at the NIT are driven by both the HOMO- 
LUMO overlap and the hydrogen bonds. Quite many in- 
termolecular contacts involving chlorine atoms have also 
been evidenced. Moreover, we have found a direct rela- 
tionship between the relative variations of the potential 
energy density at the critical point and the distance be- 
tween the atoms involved in the intermolecular bond. 

Three estimates of the charge transfer have been 
obtained from our tight-binding model, from integra- 
tion over the atomic basins defined within Bader's 
approach and from a simple model based on isolated 
molecules calculations. All three estimates are in close 
agreement with each other and give a value of about 
0.45 e" at 300 K and 0.65 e" at 40 K. Experimental 
charge density collection on TTF-CA is under progress 
in order to check the rather high value obtained at 300 K. 

All these results provide a possible picture for the 
mechanism of the NIT. The energy dispersions in VB 
and CB are mainly along a* and b* with an indirect gap 
between T and Y. Above the transition temperature, the 
VB has dominant D character at Y and the CB a pure A 
character at T. When temperature is lowered, the slow 
decrease of the cell parameters induces a slow variation 
of all transfer integrals and the band gap progressively 
vanishes. There is a point where the indirect band gap 
closes. Pockets of holes at Y and electrons at T appear 
which activate a new mechanism assisted by intraSSiSi 
and intermolecular 37,38 vibrations creating DA pairs (and 
longer . . .DADA. . . strings) for which the inversion sym- 
metry is locally broken. This picture is consistent with 
the observed pretransitionnal effects^ described as ionic 
segments in the neutral phase 3 * 3 ^*!^ sometimes called 
CT exciton strings4i The rapid increase of CT along the 
chains may then activate the hydrogen bonds via intra- 
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molecular electron-electron interactions between the va- 
lence 7r and deep a molecular orbitalsii^ This enables 
the whole crystal to topple into the ionic state at the 
transition point and is consistent with the sudden de- 
crease of the cell parameter b» 4 Once the whole crys- 
tal has switched the band gap opens and both VB and 
CB have strong D-A character leading to a macroscopic 
change of the CT. 
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